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^ . ABSTRACT 

OO ! 

The active galactic nucleus PG 1553+113 was observed by the MAGIC telescope in July 2006 during a multiwavelength campaign, 
• • , in which telescopes in the optical, X-ray, and very high energies participated. Although the MAGIC data were affected by strong 
^ atmospheric absorption (calima), they were analyzed after applying a correction. In 8.5 hours, a signal was detected with a 

significance of 5.0 cr. The integral fiux above 150 GeV was (2.6 ± 0.9) ■ 10"'' phs~^m~^ By fitting the differential energy spectrum 
KN( with a power law, a spectral index of —4.1 ± 0.3 was obtained. 

' Key words. Gamma rays: observations - Galaxies: active - BL Lacertae objects: individual: PG 1553+113 



1. Introduction 

The Major Atmospheric Gamma-ray Imaging Cherenkov 
(MAGIC) telescope, located on the Canary Island of La 
Palma at 2200 m a.s.l., is capable of extending very high 
energy observations to energies previously unreachable and 
detecting new sources at energies down to 50 GeV. 

One of these sources is the BL Lac type object 
PG 1553+113, observed for the first time in this en- 
ergy range in April and May 2005 by the MAGIC tele- 
scope and the High Energy Stereoscopic System (H.E.S.S.). 
From these observations a faint signal was measured, 
and the detection was c onfirmed by furthe r observations 
(|Aharonian et all l2006t lAlbert et all 120071) . In the sub- 



sequent years, additional VHE data were taken. From 
the combined data sets, stro ng signals were found : for 
H.E.S.S. at 10.2 cT significance (lAharonian et all [20081 ) and 
for MAGIC at 15.0 cr significance jDornerl [2008) . In addi- 
tion, the VHE measurements allowed the unknown red- 
shift of the source to be constrained. Until now the 
redshift of PG 1553+113 has been unknown, since nei- 
ther emission or absorption lines could be found, nor 
the host galaxy resolved. Several lower limits were deter- 
mined (Cara ngelo et al. 2003; Sb aruf atti et al. 2005, 20QBi 
IScaroa et aL l l2000l:TTreves et a"l.ll2007l) . With the MAGIC 
and H.E.S .S. measurements, upper limits could be de- 
termined (lAharonian et al.l 120061 : iMazin fc Goebell l2007t 
lDorneill2008l b 
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To study the spectral energy distribution of a source, 
simultaneous data from different wavelengths are manda- 



2 



J. Albert et al.: MAGIC Observations of PG 1553+113 during a MWL Campaign (RN) 



tory. Therefore, a multiwavelength (MWL) campaign was 
performed in July 2006 to observe PG 1553+113. This pa- 
per concentrates on the data taken by MAGIC during this 
MWL campaign. 

2. Observations and Data Quality 

Between April 2005 and April 2007, MAGIC observed 
PG 1553+113 for a total of 78 hours. Part of these data 
were acquired during a MWL campaign in 2006 July with 
the H.E.S.S. array of lACTs, the X-ray satellite Suzaku 
and the optical telescope KVA. Suzaku observed the source 
between 24 July, 14:26 UTC and 25 July, 19:17 UTC, and 
H.E.S.S. between 22 July and 27 July. From the KVA, data 
between 21 July and 27 July are available. 

The MAGIC telescope observed PG 1553+113 between 
14 July and 27 July for 9.5 hours at zenith distances be- 
tween 18° and 35°. The data were acquired in wobble mode, 
where the source was tracked with an offset of ±0.4° from 
the center of the camera, which enabled simultaneous mea- 
surement of the source and the background. 

One hour of data was excluded due to technical prob- 
lems. The quality of the entire data set acquired during 
the MWL campaign was affectd by calima, i.e. sand-dust 
from the Sahara in the atmosphere. For the affected data, 
the nightly values of atmospheric absorption ranged be- 
tween 5% and 40%. To account for the absorption of the 
Cherenkov light, correction factors were calculated and ap- 
phed to the data (see Sect.H]). 

3. Analysis 

The data were processed by an automated ana lysis pipeline 
(|Dorner fc Bretzll2005t [Bretz fc DorlIall2005D at the data 
center in Wiirzburg. The analysis includ es an absolute cal- 
ibration with muons (iGoebel et al. ''20051 , an absolute mis- 
pointing correction (|Riegel et al.ll2005 l. and it uses the ar- 
rival time information of the pulses of neighboring pixels 
for noise subtraction and background suppression. 

In determining the background, three OFF regions were 
used, providing a scale factor of 1/3 for the background 
measurement. 

To suppress the background, a dynamical cut in Area 
(Area=7r-WIDTH-LENGTH) versus SIZE and a cut in d 
were applied. More d etails of the cuts can be found in 
iRiegel fc Bret j ()2005[ ). a nd the aforem entioned image pa- 
rameters are described bv lHillasI ()1985D . To account for the 
steeper spectrum of PG 1553+113, the here presented anal- 
ysis applied a cut in Area that was less restrictive at lower 
energies compared to the standard cut used by the auto- 
mated analysis, which has been optimized for Crab Nebula 
data over several years. 

In generating the spectrum, the cut in Area was selected 
to ensure that more than 90% of the simulated gammas 
survived. To study the dependency of the spectral shape on 
the cut efficiency, a different cut in Area with cut efficiencies 
of between 50% and 95% for the entire energy range was 
applied. 



4. Correction of the Effect of Calima 

Calima, also known as Saharan Air Layer (SAL), is a 
layer in the atmosphere that transports sand-dust from the 



Sahara in a westerly direction over the Atlantic Ocean. The 
SAL is u sually situated between 1.5 km and 5.5 km a.s.l. 
()Dunion fc S. 2004). Since the Canary Islands are close to 
the North African coast, the MAGIC observations are prob- 
ably affected by additional light absorption in the atmo- 
sphere when calima occurs. 

From extinction measure ments of the 

Carlsberg Meridian Teles c ope ((Evans et al.l 120021 : 
ICarlsberg Meridian Telescope! |2007[ ). which are avail- 
able for each night, the loss of light due to calima was 
calculated. To correct the data for absorption by the 
atmosphere, the calibration factors were adapted for each 
night separately and the data wer e reprocessed . Mor e 
details of the method are provided bv lDorner et al.l ()2008f ). 



5. Results 

The 8.5 hours of data from PG 1553+113 provided a signa l 
with a significance of 5.0 cr according to iLi fc Mai (|T983l ). 
The d"^ distributions for the ON- and normalized OFF- 
source measurement are shown in Fig. [TJ 
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Fig. 1. Distributions of d"^ for ON-source events (black 
crosses) and normalized OFF-source events (gray shaded 
area) from 8.5 hours of data. The dashed line represents 
the cut applied in d"^ . 



The differential spectrum measured by MAGIC is shown 
in Fig.[2l In this plot, a set of spectra obtained with different 
cut efficiencies and different simulated spectra is indicated 
by a gray band. The data points of the spectrum are given 
in Table [T] including the statistical er rors. The sys tematic 
errors of the analysis are discussed in I Albert et al.l (^008). 

Fitting the differential spectrum with a power law yields 
a flux of (1.4 + 0.3) • 10-s ph TeV-^s^^m-^ at 200 GeV and 
a spectral index of —4.1 + 0.3. This result is consistent with 
the data (fit probability 45%) and in good agreeme nt with 
previous measurements in 2005 and 2006 ( Aharonia n et al.l 
l2006l : I Albert et al.l[2007i l. Within the errors, the simuhane- 
ous measurements of the H.E.S.S. telescopes in the energy 
range above 225 GeV is in agreement with the MAGIC re- 
sults, although the fit of the differential spectrum yields a 
spectral index of -5.0 ± 0.7 (Aharonian et al. 2008). The 
integral flux above 150 GeV obtained from this analysis is 
(2.6 ±0.9) • 10"^ phs-im-2 
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Fig. 2. Differential energy spectrum of PG 1553+113. The 
fiorizontal error bars indicate the width of the energy bins. 
The vertical error bars illustrate the statistical errors. The 
gray band corresponds to a set of spectra obtained with 
different cut efficiencies and simulated spectra. 
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Fig. 3. Night-by-night light curve for the integral flux above 
150 GeV in phs~^m~^ (upper part) and the flux in the 
optical R-band (lower part) between July 15*'' and 27*'' 
2006. For the first nights of the MAGIC observations, no 
data by the optical telescope KVA were obtained. 



Table 1. Flux for the spectral points in Fig. [5] including 
the statistical errors. 



E 
[GeV] 


F 

[phTeV'^s^^m-^] 


Statistical Error 
[pliTeV-^s-^m-2] 


91 


3.36-10-'' 


8.11-10-" 


139 


7.37-10-^ 


2.68-10-'^ 


211 


1.67-10-^ 


6.04-10-^^ 


320 


9.61-10-* 


9.61-10-* 


486 


6. 18-10-* 


4.00-10-* 



Table 2. Date, start time, duration, and flux above 
150 GeV for the MAGIC observations carried out between 
July 15*'' and 27*'' 2006. For one day, the flux level could not 
be determined because no signal for energies above 150 GeV 
was detectable in this data set. 



Observation [UTCJ 
Start Time (Duration) 



F (>150GeV) 
[phs-^m-^] 



15.7.2006 21:45 (1.4 h) 
16.7.2006 21:45 (1.37 h) 
17.7.2006 22:15 (0.85 h) 
21.7.2006 22:13 (0.62 h) 
23.7.2006 21:38 (0.88 h) 
24.7.2006 21:50 (1.09 h) 
25.7.2006 21:38 (1.31 h) 
26.7.2006 21:37 (0.92 h) 



(1.55 + 1.46) - 10" 
(4.79 + 2.31) - 10" 
(3.47 + 2.14) - 10- 
(2.99 + 2.04) - 10" 

(3.34 + 2.38) - 10" 
(1.50 + 2.18) - 10" 
(1.27 + 1.87) - 10' 



To check whether a flare occurred during the eight 
nights of observation, the flux above 150 GeV was calcu- 
lated on a night-by-night basis (see Table [2|) . The corre- 
sponding light curve is shown in the upper part of Fig. [3l 
Since the data sets for single days are of durations shorter 
than one and a half hours (0.62h- 1.4h), the data points 
represent only signals of a significance level between 0.7 <t 
and 2.1(7. Consequently, no strong conclusions about the 
night-to-night variability in the flux can be drawn. Within 
the errors, the measurement is consistent with a constant 
flux (fit probability 82%). 



Contemproraneously with the MAGIC observations, the 
optical telescope KVA acquired data in the R-band. For the 
first nights, no data was available. The flux for additional 
nights is shown in the lower part of Fig. [31 



6. Conclusion and Outlook 

PG 1553+113 was observed in July 2006 as part of a MWL 
campaign with the MAGIC telescope. After correcting the 
data for the effect of calima, the analysis detected a gamma- 
ray signal of 5.0 standard deviations. Within the statistical 
errors, the differential energy spectrum is compatible with 
those derived by previous measurements including the one 
observed by H.E.S.S. in this campaign. The inter-night light 
curve shows no significant variability. 

The measured fiux and reconstructed spectrum will be 
used in studies of the spectral energy distribution, which 
will include other da ta acquired during the MWL campaign 
()Reimer et al.ll2008[ ). 
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